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We develop a technique, denoted as the finite radius approximation (FRA), that uses a two- 
dimensional version of the Shannon-Nyquist sampling theorem to determine transverse densities 
and their uncertainties from experimental quantities. Uncertainties arising from experimental un- 
certainties on the form factors and lack of measured data at high Q 2 are treated. A key feature 
of the FRA is that a form factor measured at a given value of Q 2 is related to a definite region 
in coordinate space. An exact relation between the FRA and the use of a Bessel series is derived. 
The proton Dirac form factor is well enough known such that the transverse charge density is very 
accurately known except for transverse separations b less than about 0.1 fm. The Pauli form factor 
is well known to Q 2 of about 10 GeV 2 , and this allows a reasonable, but improvable, determination 
of the anomalous magnetic moment density. 
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I. INTRODUCTION 

A truly impressive level of experimental technique, 
effort and ingenuity has been applied to measuring 
the electromagnetic form factors of the proton, neu- 
tron (nucleon) and pion drill- These quantities are 
probability amplitudes that a given hadron can ab- 
sorb a specific amount of momentum and remain in 
the ground state, and therefore should supply infor- 
mation about charge and magnetization spatial den- 
sities. 

The text-book interpretation of these form factors 
is that their Fourier transforms are measurements of 
the charge and magnetization densities. This inter- 
pretation is deeply buried in the thinking of nuclear 
physicists and continues to guide intuition, as it has 
since the days of the Nobel prize-winning work of 
HofstadterQ. Nevertheless, the relativistic motion of 
the constituents of the system causes the text-book 
interpretation to be incorrect Q- The difficulty is 
that in electron-proton scattering the initial and final 
nucleon states have different momenta and therefore 
different wave functions. In general, these different 
states are related by a boost operator that depends 
on the full complexity of QCD. The use of transverse 
densities [H [l(| avoids this difficulty by working in 



the infinite momentum frame and taking the space- 
like momentum transfer to be in the direction trans- 
verse to that of the infinite momentum. In this case, 
the different momenta of the initial and final nucleon 
states are accommodated by using two-dimensional 
Fourier transforms. The transverse charge and mag- 
netization densities are constructed from density op- 
erators that are the absolute square of quark-field 
operators, so they are correctly defined as densities. 

This paper is concerned with extracting the spa- 
tial information by developing and using a theoretical 
technique that is model-independent and also pro- 
vides a practical way of dealing with both experi- 
mental uncertainties and the lack of information on 
unmeasured regions, with minimal assumptions. In 
the subsequent text we plan to show how to construct 
bands of transverse densities that are consistent with 
available experimental knowledge and also take into 
account the possible effects of data taken at momen- 
tum transfer Q 2 higher than available in the present 
data set. This allows one to consider the possible 
impact of future experiments. 

But there also is a more general context, with 
the high current interest in mapping the three- 
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dimensional structure of the nucleon [11| . Therefore 
we also aim to provide a technique that can be eas- 
ily extended determining the spatial aspects of other 
quantities related to transverse momentum distribu- 
tions and generalized parton distributions. 

Next we present an overview of the remainder of 
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this paper. Sect. |TT] concerns the following situation. 
Suppose a form factor F{Q 2 ) and transverse density 
p(b) are related by a two-dimensional Fourier trans- 
form, and that p(b) is localized, pip) = for b greater 
than some finite distance. The function p(b) is band 
limited and can be written as a discrete Fourier series 
involving F(Q 2 ). This result, known as the Nyquist- 
Shannon [12J sampling theorem, enables us to asso- 
ciate the density at a given range of values of b with a 
discrete value of the momentum transfer, see Eq. ([3]) 
below (which we denote as the finite radius approxi- 
mation FRA). The equivalence between the FRA and 
the Bessel series expansion technique is also estab- 
lished. A general version of the FRA, applicable to 
other observable quantities, is also presented. 

Sect. IIIII is concerned with exploring the the valid- 
ity and utility (which depends on the number of terms 
needed in the discrete Fourier series) of the FRA us- 
ing examples in which the form factor is given by a 
monopole (M) or dipole (D) form. Sect. [V] is con- 
cerned with the reality that the proton electromag- 
netic form factors are not known as analytic func- 
tions. Instead, form factors Ge,m, F\2 (with uncer- 
tainties) measured at discrete values of Q 2 up to a 
finite maximum value Q^ax are known. This means 
that p is known only within some uncertainties, and 
a technique to determine the uncertainties in p must 
be developed. This is accomplished by using the val- 
ues of Fi ± dFi in the FRA. Estimates of the effects 
of incompleteness, arising from contributions in the 
unmeasured region, Q 2 > Qm ax > are a ^ so provided. 
The paper is concluded with a brief summary. 

II. GENERAL CONSIDERATIONS 

Intuitively, we expect particles to be localized. 
That is, we expect densities associated with the par- 
ticle to be well approximated by functions that are 
zero outside some maximum radius. This assump- 
tion, called the finite radius approximation (FRA), 
greatly simplifies the relationship between form fac- 
tors and their associated densities. 

Let p(b) be a two-dimensional transverse density 
function (we later take this to be charge or magneti- 
zation density) and let F(Q 2 ) be the associated form 
factor. The transverse density is given by [H, [l3[ 

P(b) = (2^5) / dV- 4q b F(Q 2 = q 2 ) 

= h J Q d QMQb)F(Q 2 ), (1) 

with the azimuthal symmetry of p obtained from the 
Lorentz invariant form of F in the space-like region 



with q + = 0. If one knows F(Q 2 ) exactly for all 
values of Q 2 , the transverse density is known im- 
mediately. However, one only knows F(Q 2 ) within 
experimental uncertainties for a finite range of Q 2 . 
This means that p is known only within some uncer- 
tainties, and it is necessary to develop a technique to 
determine the uncertainties in p. 

We proceed by assuming that p(b) ps for b > R, 
where R is a finite distance. Since the functions p, F 
are Fourier transforms, F is band-limited. We pro- 
ceed in the spirit of the Nyquist-Shannon sampling 
theorem. The function p can be expanded as 

/>(&) = $>„ Jo (X„-), (2) 

n— 1 

where X n is the n-th zero of Jo, and c„ is given ap- 
proximately by the formula 

1 2 

C "~ C " = 2^R 2 J 1 (X n ) 2F{Q2nh (3) 

with 

Qn = ^. (4) 

The above equation Eq. ([3]), which is the two- 
dimensional version of [12(, is the central formal re- 
sult of this paper. Using this in Eq. ([2]) yields the 
following expression for p(b): 

1 00 h 

n— 1 

The result Eq. ([5]) is the central phenomenological re- 
sult because it tells us that measuring a form factor at 
Q 2 n provides information about the density mainly at 
values of b < RjX n . This is because Bessel functions 
are of the order of unity only for values of arguments 
less than that of its first zero. 



A. Equivalence with the Bessel Series 

Replacing c„ by c„ would be exact if the assump- 
tion p(b > R) = is exactly true. This condition 
is clearly approximately true, so we expect a near 
equality between c n and c n . In fact, it turns out that 
the approximation is amazingly accurate as we now 
demonstrate. Numerical examples are provided in 
subsequent sections. The exact values of c n are ob- 
tained from the orthogonality of the cylindrical Bessel 
functions as 

c. = ^^f b P (b) M x n t )db . (6) 
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The use of this in Eq. (fTJ) followed by integration over 
b can be done using a standard identity to yield 



Xn r qF(q 2 )J (qR) 

^E?J!{X n ) J {^) 2 -q 2 



dq. (7) 



We may use a dispersion relation for the form fac- 
tor [l4j to establish the connection between c„ and 
c„. First recall that, for Q 2 > 0, 



F(Q 2 



1 



dt 



ImF{-t) 



(8) 



T J Ami 1 + Q 

and that using this expression in Eq. ([l} yields 
1 f°° 

P( b ) = 7TS dt K (Vib)ImF(-t). (9) 

^ J Ami 

Proceed by using the above in Eq. © and then inte- 
grate over b using: 



bdbK (Vtb)J (b 



X n 

~r 



1 



+ t 



1 + X n ,h{X n )K Q {VtR) 



(10) 



Then, 



R*Jl{X n )W 



(II 



1 ///-. 



ImF(-t) 



l + X n ,h{X n )K {VtR) 



(11) 



Using only the first term within the brackets along 
with Eq. ((SJ allows one to identify the integral over 
t as nF(Q 2 l ). Thus (using Eq. ((3])) one arrives at the 
result that c n — c n plus a correction term, suppressed 
by a modified Bessel function evaluated at a large 
argument. For example, a significant contribution to 
the Im F comes from the region t ~ m 2 = 0.5 GeV 2 , 
and using R = 3.3 fm (see Sec. then K (m p R) = 
10~ 6 . The net result is that 



J{(XnW J Ami t+ljp 



A reasonable estimate is that 

c n c 5 



X n Ji{X n )W 



(12) 



(13) 



The condition that S n be small is that R be chosen 
to be large enough. We ensure that this condition is 
well-satisfied for all of our examples and applications. 



B. Preliminary Evaluations 

It is worthwhile to perform some preliminary anal- 
ysis of the expression Eq. ((5]). For x 3> 1, Jo(x) is 
well approximated [la] by 



. , /2 , 7T. 

Jq(x) » W — cos(x- -) 

V 7TX 4 



(14) 



so that the n'th zero of Jo, X n , is given approximately 

by 



X n w (n+ -)tt , 



(15) 



and 



•/l(^n) — ~«^o(^n) 

ps (-l)"2 1 / 2 vr- 1 ((n- 



4 yy 



(16) 



It follows that for large n, the terms in the series 
Eq. ([5]) for p(b) are of the form: 



7T 



2i? 2 



^)F(Q 2 )J (X„A) 



at b = 0. So for the series to converge everywhere, 
namely at b=0, we need F to fall faster than Q~ 2 
for large Q. The oscillations of the cylindrical Bessel 
functions hastens the convergence for non-zero values 
of b. 

Given this convergence, the function p(b) can be 
approximated by using a finite number of terms in 
the series Eq. Because Q 2 = (X n /R) 2 serves as 
the Q 2 in the argument of F, cutting off the series 
at N terms is equivalent to taking F(Q 2 ) = for 
Q 2 > (X N /R) 2 . 

If the assumption that p(b) — for b > R holds for 
a given value of R, then it also holds for larger values 
of R. We can see from Eq. ([5]) that increasing R 
increases the frequency with which F{Q 2 ) is sampled 
and therefore decreases the range that is sampled. As 
a consequence, an increase in R demands an increase 
in the number of terms in the approximation for p. 

A quick result following from the fact that p is the 
Fourier transform of F is that the mean-square-radius 
(b 2 ) is given by 



(b 



2\ — 



J d 2 b b 2 p(b) = -4 



dlogF 



dQ 2 



Q 2 =0 



(17) 



In this paper, we choose R ~ 5-\/|(& 2 )| in determin- 
ing the number of terms in our expansion. Numerical 
studies of the form factors considered in preparing 



4 



this paper have shown that this value of R is suffi- 
ciently large so that perturbations to this value lead 
to the same density functions and that R 2 p(R) is al- 
ways small enough so that the difference between c n 
and c n is minute. 



C. Other transverse densities 

We believe that the techniques used in this pa- 
per can be exploited to image other quantities that 
depend on transverse position. Suppose there is a 
transverse quantity p^ (b) that is a two-dimensional 
Fourier transform of an experimental observable 
FW(Q 2 ) such that 

P W (b) = i- J QdQJ x {Qb)F {x) (Q 2 ). (18) 

An example, discussed in detail in Sect. IV CI is the 
magnetization density p m of the anomalous magnetic 
moment. We expect that the index (A) is associated 
with a given number of units of orbital angular mo- 
mentum. Extracting p( x >(b) is facilitated by using 
the expansion 

oo , 

p W(b) = Y,CnxJx(X x , n j l ), (19) 

n— 1 

where X\ >n is the n'th zero of the Bessel function of 
order A. Then the sampling theorem leads immedi- 
ately to the result. 

C «.A « C„, A = 2 fW(Q| ) ; (2Q) 

R z j\ + i(Xx, n y 

n _ X\^ n 



The difference between c n \ and c nt \ can be shown to 
be very small by using the arguments of Sect. Ill Al 
The result Eq. (|2T))) can be used to relate accessible 
kinematic ranges with transverse regions. 

III. EXAMPLES 

To demonstrate our method and explore its limita- 
tions, we now analyze two models of the form factor. 
For the first model, let the form factor be given by 
the monopole form 

F M {Q 2 ) = (21) 

where f2 = 0.77 GeV. This form factor is taken as 
a caricature of the pion electromagnetic form factor. 



Then the associated charge density is obtained from 
Eq. (TTJ: 

PM(b) = ^-n 2 K {nb). (22) 

This function diverges as log (1/6) for small values of b 
and so provides a severe test of the method. With the 
stated value of Q we find (b 2 ) M = 4/fl 2 = 0.26 fm 2 , 
and thus take R = 5\/\(b 2 )\ = 2.56 fm. We then 
find the fractional difference between c„ and c n of 
Eq. (TT2"|) is less than 5xl0 -4 for small values of n, 
and the magnitude decreases rapidly as n increases. 
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FIG. 1: (Color Online) Plot of p M (blue,solid), 10 term 
approximation (red, long dash) and 20 term approxima- 
tion (green, medium dash) and 50 term (brown, short 
dash). 

We compare to pu to its approximation as an ex- 
pansion in N terms, with N = 10, 20, 50 in Fig. [TJ 
We see that our approximations differ from the exact 
result, but the difference decreases with increasing 
value of N. The 50 term approximation works rea- 
sonably well for all value of b for which the density 
differs appreciably from 0. Unfortunately the 10,20 
and 50 term approximations would require measure- 
ments at Q 2 = 6,23 and 144 GeV 2 . Only the first 
value seems presently achievable. 

We now examine the dipole form factor given by 

^=(TTfjs (23) 

where A 2 = 0.71 GeV 2 . This value is suggested by its 
historically close relationship with the proton electro- 
magnetic form factors. The dipole transverse charge 
density is obtained by from Eq. (JTJ) to be 

p D {b) = -Ua^&A). (24) 

47T 

This form factor falls more rapidly with increasing 
Q 2 than does Fm , and also corresponds to the larger 
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physical extent of the proton as compared to the 
pion. Furthermore, pp is not singular at the origin 
(~ 1 — 0.058(6A) 2 ). Thus there are several reasons 
to expect to find better convergence properties, and 
therefore a more accurate representation of the trans- 
verse density for the proton. With this value of A, 
(b 2 ) = 8/A 2 = 0.439 fm 2 , and R = 3.31 fm. Once 
again the fractional difference of Eq. (|12p is truly tiny 
for all values of n: the fractional differences are less 
than about 10 -5 for all values of n that correspond 
to non-zero c„. We plot po and its approximations 
in Fig. H 




b [fm] 



FIG. 2: (Color online) Plot of po (solid), 5 term approx- 
imation (red, long dash), 10 term approximation (green, 
medium dash) and 15 term approximation (brown, short 
dash) . 

We can see how the approximations converge to 
the exact po- Even the 10 term approximation is 
reasonably good and the 15 term approximation is 
extremely accurate except for b < 0.1 fm. 

Another way of looking at convergence properties 
is to examine properties of the transverse density. We 
display upper limit dependence of (6 2 )m,d for both 
the monopole and dipole form factors. We compute 
these matrix elements for a range of values of b from 
to 1.5 fm. This covers the region up to where p 
is about 0.1% of its central value. The results are 
shown in Table [U Despite the relatively poor conver- 
gence obtained for the monopole form factor (Fig. [1} , 
reasonable convergence for the expectation value is 
obtained. However, the convergence is much better 
for the dipole form factors. The 5,10 and 15 term 
approximations correspond to values of Q 2 = 0.9,4 
and 9 GeV 2 . These values and even higher have al- 
ready been achieved experimentally. Thus we rea- 
sonably expect that the proton transverse density is 
now known. Indeed, this has already been suggested 
[9]. However, now we can answer the question: "How 
well is the proton transverse charge density known?" . 



TABLE I: Upper limit, N, dependence of (b 2 )M,D com- 
puted for values of b from to 1.5 fm. 



N 




(b 2 )o [fmf 


5 


0.259 


0.313 


10 


0.362 


0.320 


15 


0.368 


0.319 


oo 


0.367 


0.319 



IV. EXTRACTION OF PROTON FORM 
FACTORS AND UNCERTAINTIES 

The transverse densities we seek are given in terms 
of the Dirac Fi and Pauli F% form factors, which are 
expressed in terms of the Sachs electromagnetic form 
factors Ge and Gm as 



Fi{Q 2 



where r 



G t 



tG 



M 



1 + T 



,F 2 (Q 2 



G 



M 



G h 



(25) 



Elastic electron-proton scattering has been mea- 
sured up to Q 2 of about 30 GeV 2 , with the separa- 
tion of both G e and Gm extracted using a variety of 
techniques up to 10 GeV 2 . There are two sources of 
uncertainty in the extraction of the transverse den- 
sities. Experimental uncertainties from the measure- 
ments of Ge and Gm yield uncertainty in the ex- 
tracted densities, and incompleteness error arise from 
the lack of form factor measurements at very high Q 2 
(above 30 GeV 2 ). In this section, we perform extrac- 
tions of the transverse density and evaluate the the 
effects that two kinds of uncertainties on the densi- 
ties. 

The form factors Ge and Gm have been extracted 
from a global analysis of the world's cross section 
and polarization data, including corrections for two- 
photon exchange corrections from Ref. [l6j . The 
analysis is largely identical to that that of Ref. UJl . 
although additional high Q 2 form factor results [l8| 
have been included. In addition, the slopes of Ge 
and Gm at Q 2 =0 were constrained in the global fit 
based on a dedicated analysis of the low Q 2 data. In 
the global fit, the large body of high Q 2 data, espe- 
cially for Gm, can constrain the fit well enough that 
the low Q 2 behavior is not primarily constrained by 
the low Q 2 data. Constraining the slope based on an 
analysis of only the low Q 2 data keeps the global fit 
from doing a poor job at low Q 2 simply to make a 
slight improvement in the high Q 2 data. In writing 
Ge{Q 2 ) = 1 — Q 2 Re/6, the value of Re was con- 
strained to be 0.878 fm and Rm was constrained to 
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be 0.860 fin. This is important in the extraction of 
the large scale structure of the density. The fit is of 
the following form: 



G M (Q 2 ) =Hj, 
G e (Q 2 ) = 



I+Pgt+Piqt 2 +P14T 3 



1 + p 2 T + p 4 T 2 + p 8 T 3 + p 12 T 4 + p 16 T 5 

l + q(jT + q 10 T 2 + qnT 3 



1 + q 2 T + q^T 2 



qi2T 4 + qi 6 r 5 



(26) 



where the fitting constants pi, ..pi6, Q2, S16 are 
given in Table [TTI and we use \i p = 2.792782. 



TABLE II: Fit parameters for Gm (pi), G E (qi[ 



i 


Pi 


9i 


2 


9.70703681 


14.5187212 


4 


3.7357 x 10~ 4 


40.88333 


6 


-1.43573 


2.90966 


8 


6.0 x 10" 8 


99.999998 


10 


1.19052066 


-1.11542229 


12 


9.9527277 


4.579 x 10~ 5 


14 


2.5455841 x 10" 1 


3.866171 x 10~ 2 


16 


12.7977739 


10.3580447 



We also need a reliable estimate of the experimen- 
tal uncertainties in the form factors, in order to de- 
termine the uncertainty in the extracted coefficients 
c„ . In the global analysis, there are two sources that 
can contribute to the uncertainties in Ge and Gm- 
the uncertainty on each individual cross section or 
polarization ratio, and the normalization uncertainty 
associated with each cross section data set. The nor- 
malization factors are allowed to vary in the fit, as 
was the case in Ref. 17]. To estimate the uncertainty 
in the fitted normalization factors, we take the nor- 
malization factor from a single data set and vary it 
around its best fit value (while allowing all other pa- 
rameters to vary) to map out the change in the \ 2 of 
the fit as a function of the normalization factor. This 
yields uncertainties between 0.2% and 2.5% (typically 
0.6%-l%), compared to the initally quoted uncer- 
tainties of 1.5% to 5%, for the data before the nor- 
malization has been constrained by the fit. However, 
by assuming that all uncertainties are entirely un- 
corrected or pure normalization factors, we neglect 
the possibility there may be some angle-dependent 
or Q 2_ dependent correction that could bias the de- 
termination of the relative normalization coefficients. 
Thus, we assume that the final uncertainty on each 
normalization factor is at least 0.5%, even if the re- 
sult of the x 2 analysis yields a smaller result. 



Having the uncorrelated uncertainties for each data 
point and the constrained normalization uncertain- 
ties, we then extract the uncertainties for Ge and 
Gm- For the uncorrelated uncertainties, we randomly 
shift each cross section and polarization ratio mea- 
surement within its uncertainties, and then redo the 
fit for Ge and Gm- We repeat this 1000 times, and 
look at the range of values for several Q 2 values (55 
Q 2 values between 0.007 and 31.2 GeV 2 ). This yields 
our uncorrelated uncertainty at each of the Q 2 points. 
To obtain the impact of the normalization uncertain- 
ties, we repeat this procedure, varying the normal- 
ization of each cross section data set according to its 
uncertainty, and determine the range of Ge, Gm val- 
ues for the same set of Q 2 points. In this procedure, 
the uncertainty obtained depends on the fit function 
used, and a functional form with insufficient flexibil- 
ity will yield significant smoothing of the results and 
thus unrealistically small uncertainties. We scale up 
our uncertainties by a factor of two, which yields good 
agreement with best direct measurements of the form 
factors and uncertainties. 

As mentioned above, we use the electric and mag- 
netic radii extracted from just the low Q 2 data in as 
a constraint to the global fit, which can yield unreal- 
isically small uncertainties for below Q 2 =0.2 GeV 2 , 
especially for Gm , where the very low Q 2 data is ex- 
tremely limited. Thus, for these low Q 2 values, we 
calculate the uncertainty at each Q 2 corresponding 
to the uncertainty in the extracted radius, assuming 
the linear expansion. We take this larger uncertainty, 
rather than the result from the fit, until the uncer- 
tainties from direct extractions of the form factors are 
of comparable size, at which point we take the direct 
extraction of the uncertainty. For Q 2 > 10 GeV 2 , 
there are no direct extractions of Ge, and thus we 
again have to be sure that we do not underestimate 
the uncertainties. The global fit yields Ge/Gd ~ 
at high Q 2 , but it is difficult to tell if Ge becomes 
zero, or if Ge / Gm continues its linear decrease with 
Q 2 [ll]. Thus, for Q 2 > 10 GeV 2 , we set the uncer- 
tainty to be the difference between the best fit, which 
yields Ge ~ and the fit where the linear falloff in 
Ge/Gm continues, with Ge changing sign and then 
increasing in absolute value. 

We then use the fit and uncertainties for Ge and 
Gm to extract F\ and F 2 , treating the uncertainties 
in Ge and Gm as uncorrelated, yielding: 



1 



:) 2 (dG M ) 2 



(dF 1 y = {—y(dG E y + (- 



(dF 2 ) 2 = (-—) 2 (dG E ) 2 + {-—) 2 {dG M f 



(27) 



(28) 
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FIG. 3: (Color online) The electromagnetic form factors 
Fi((5 2 ) and F 2 {Q 2 ) and their error bands, scaled by a 
factor of Q 4 . 



While the Rosenbluth extractions yield a strong anti- 
correlation between the uncertainties onGg and Gm, 
the polarization ratio yields a correlated uncertainty; 
in the global fit, the combined result is fairly well 
approximated by entirely uncorrelated uncertainties. 
Figure [3] shows the extracted values of F\ and F 2 
along with their uncertainties. Because the elastic 
cross section is dominated by the contribution from 
Gm at large Q 2 , the fractional uncertainties on Ge 
are much larger, and the uncertainty on domi- 
nates the uncertainty on both F± and F2, even though 
its contribution to Fi is supressed by a factor of r rel- 
ative to the Gm contribution. 

We note that for Q 2 < 0.5, the uncertainty coming 
from cross section normalizations can be the larger 
contribution to the total uncertainty (and it's dom- 
inant for Ge below 0.1 GeV 2 ). While the normal- 
ization uncertainty in the cross sections won't give 
a normalization style uncertainty on Ge, the nor- 
malization of a given experiment will tend to have 
a correlated effect on all of the extractions within 
the Q 2 covered by the experiment. This effect is ac- 
counted for by using the procedure discussed below 
in Sect. [O] 



V. EXTRACTION OF REALISTIC PROTON 
TRANSVERSE DENSITIES 

The principle aim of this paper is to use data 
observed in experiments to obtain the charge and 
magnetization densities. Recall that the transverse 
charge density p c h is given by 

Pch(b) = QdQMQtyF^Q 2 ). (29) 

The two-dimensional Fourier transform of F%, p 2 is 
similarly given by 

p 2 (b) = ^f QdQJ (Qb)F 2 (Q 2 ). (30) 

However the true magnetization density, obtained by 
computing the expectation value of the transverse po- 
sition operator with the electromagnetic current op- 
erator is given [To} by 

Pm{b) = -b-^p 2 (b) 

= h J Q 2d Q J i(Qb)F 2 (Q 2 ). (31) 

This quantity is the density related to the anoma- 
lous magnetic moment. We begin by extracting p c h,2- 
The starting point is to use the above expressions 
along with the experimentally determined F\ >2 ob- 
tained from the fits of Sect. IIVI But extracting re- 
alistic transverse densities requires that a determina- 
tion of the uncertainties in the results. There are two 
sources of uncertainty. Experimental data have un- 
certainties in the region where they are measured, 
and no direct information is available above some 
maximum value of Q 2 = Q^^, where there are no 
measurements. The experimental uncertainties lead 
directly to uncertainties in the c„ via Eq. ([3]) , and can 
be taken into account without further ado. However, 
uncertainty must arise because of lack of knowledge 
of form factors for Q 2 > Q 2 lax , and these need to 
be estimated. This error is called the incompleteness 
error. 



A. Impact of Experimental Uncertainties on 
the Extracted Transverse Densities 

We first treat the experimental uncertainties. We 
only use the series Eq. (0 for values of Q 2 n for which 
form factors have been extracted. The magnetic form 
factor Gm is well measured up to Q 2 = 31 GeV 2 , but 
Ge is only known up to ~10 GeV 2 . Based on the 
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estimated uncertainties on Ge above 10 GeV 2 , we 
find that while F± is relatively well measured up to 

30 GeV 2 , the uncertainties on F 2 grow rapidly above 
10 GeV 2 , reaching 25% by 13 GeV 2 . These upper 
limits on Q 2 are related to limits on the summation 
index n (of Eq. (|S"])) through Eq. Q which requires 
values of Ri. Taking (b 2 ) given by p c h,2 from the fits 
presented above, we use Eq. (fTT)) to obtain Ri = 3.29 
fm and R2 — 3.62 fm for F\ 2- This corresponds to 
upper limits N on the sum over n n = 30, Q 2 — 

31 GeV 2 for p ch (b), but only up to n = 20, Q| = 11 
GeV 2 for p 2 (b). 

p c h(b) [fm 
2.0, 




b [fm] 



0.0 0.2 0.4 0.6 0. 



1.0 1.2 1.4 



b [fm] 



FIG. 4: (Color Online) The transverse densities p c h,p2 
(blue, solid) of the parameterizations and their approx- 
imates to 10 terms (red,long dash), 20 terms (green, 
medium dash), 30 terms (brown, short dash) and using 
the parameterization of Eq. (|26[) . The approximations 
converge as the number of terms increases. 



A exp p(b) [fm" 2 ] 
0.030 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



FIG. 5: (Color online) Uncertainties in transverse den- 
sities A eX pp c h (solid, blue) and A e x P p2(b) (dashed, red) 
due to experimental uncertainties on Fi, F%. 



:>,() 



obtain 

&expPch{b) = ^2 

71=1 

30 



dp c h(b) 



dFi 



Hi 



1 n— 1 



Jo{X n 



Ri 



dFiA 2 ), 



(32) 



211 



A exp p 2 (b) = 22 

71=1 

1 20 



d P2 (b) 



dF 2 



dF 2 A 2 ) 

tl2 



irR 2 

z n=l 



.h{X n — ) 
Hi 



n-2 



(33) 



Note that the errors are added linearly. This means 
that we are taking the worst case possible by assum- 
ing a full correlation. These uncertainties in densities 
are plotted in Fig. [S] They are about 1.5% of the 
transverse density at b = and decrease (in absolute 
value) at increasing distances. The fractional uncer- 
tainty is small (below 10%) until b 1 fm, where the 
density is only a few percent of the peak density. 



The transverse densities p c h,2 are plotted as the 
solid curves in Fig. 01 The densities peak at b = and 
that the transverse density p 2 has a slightly broader 
spatial extent than that of p c h- 

The next step is to extract c n from the fit to 
the form factor using Eq. ([3]). The uncertainty on 
Fx, 2 (Q 2 ) directly yields an uncertainty on c„, and 
thus its contribution to p(b) (Eq. ([2])). Assuming the 
errors from each c n extraction add constructively, we 



B. Incompleteness Error 

We next study the uncertainties in the transverse 
density caused by lack of experimental knowledge at 
large values of Q 2 . The first step is to understand 
the meaning of the truncations made in Eq. (|29|) and 
Eq. (|3"Tj) . Plots of these approximations are given 
in Fig. |U We see that for, p 2 , one achieves agree- 
ment with the parameterization for values of N as 
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low as 20, with the largest disagreement at b — 0. For 
Pch(b — 0), the difference between the result from the 
parameterization and the N — 30 approximation is 
-2%, while for p2(b = 0), the N = 20 approximation 
is only 1% below the full result. Even though fewer 
terms are included in the approximation for F%, the 
agreement is comparable, due to the more rapid fall- 
off of F2 with increasing values of Q 2 . 

A mcP (b) [fm" 2 ] 
0.04 



Pch(b) [fm 
2.0 




0.0 0.2 0.4 0.6 0.! 



1.0 1.2 1.4 



b [fm] 



FIG. 6: (Color online) Incompleteness error. The abso- 
lute error in p c h (solid, blue) and pi (dashed, red). 

Given this information, we can state our procedure. 
Our basic transverse densities are obtained by using 
the parameterization Eq. (|2"o| to evaluate the expres- 
sions of Eq. (p5|).Eq. (1301) . and Eq. (JST). However, we 
are justified in using this parameterization for values 
of Q 2 corresponding to N = 30, (20) for F 1:(2) . We 
assume a maximum error by taking the uncertainty in 
the form factor to be ± the value given by the param- 
eterization. Therefore the estimated incompleteness 
uncertainty is given by the expression 



\nc(b) 



E 

N+l 



CnMXJR^FiQl 



(34) 



as a function of b, with i — 1,2. The results are 
shown in Fig. [5] It is necessary to realize that us- 
ing this expression for the incompleteness error over- 
estimates the error because using this expression is 
equivalent to assuming that the form factor vanishes 
for Q 2 > Q 2 N in Eq. ([3]). But the form factor can not 
suddenly drop to 0. Fig. [3] shows a fractional error 
bar for Ft (31 GeV 2 ) which is only about 0.2, a frac- 
tional error bar at 13 GeV 2 which is only about 0.3 of 
the form factor F 2 . Thus using Eq. (IMl) amounts to 
making an overestimate. To be conservative, we ob- 
tain the total uncertainty by adding the contributions 
of Eq. (|32l) (or Eq. (|3"3"|l ) to the estimated incomplete- 
ness uncertainty given by Eq. (|34[) . 

We now have working expressions for the trans- 
verse densities p c h,2, and their respective uncertain- 




b [fm] 



*r b [fm] 



0.2 0.4 0.6 



1.0 1.2 1.4 



b [fm] 



FIG. 7: (Color online) p c h (solid, blue) (with error bands 
(short dashed, red)). 



ties. We start with the basic term for p c h,2, obtained 
by using the parameterization Eq. (|26[) to evaluate 
the expressions of Eq. ([29]), Eq. fl30j), and Eq. (|3"Tj) . 
then add the two separate errors Ai nc ^ eX p to get a to- 
tal error A = A inc + A exp for p c h . A band is formed 
by considering the region between the basic plus or 
minus the appropriate A for the two densities. 

The transverse densities p c h,2(b) are plotted with 
their error bands in Fig. [7] and Fig. [S] The errors 
are very small except for values of b less than about 
0.1 fm. The results in this figure are the central nu- 
merical findings of this paper. The transverse densi- 
ties are known very well indeed. The spatial extent of 
p2 is broader than that of p c h as previously observed 
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Pm(t>) 

1.2 




FIG. 9: (Color online) The true magnetization density 
p m - The uncertainties are numerically negligible. 



vanish at b = i?2- Then the FRA gives the result: 

Pm = J2 2 (Xl,n)bQl,nF 2 (Ql n ) Jl(Ql,„&), 



Ql 



i?2 



(35) 



Once again we include the effects of the experimen- 
tal error and the incompleteness error. This latter 
error is larger in this case than for p 2 because of the 
explicit factor of X\^ n . The result for p m and its er- 
ror bands are plotted in Fig. |H1 This quantity has 
a broader spatial extent than p 2 , possibly resulting 
from the importance of the pion cloud in causing the 
anomalous magnetic moment. The uncertainties on 
this quantity are greater than for the other densities. 
Future measurements extending knowledge of F 2 to 
higher values of Q 2 would reduce these higher uncer- 
tainties. 



FIG. 8: (Color online) p2, with error bands 

VI. SUMMARY 



[l9j . Note that the realistic transverse densities dif- 
fer substantially from the dipole result of Eq. (IMf . 
shown in Fig. 2. 

C. Extraction of p m (b) 

We now turn to the true transverse anomalous 
magnetic density of Eq. (j3"Pj) . defined by taking the 
matrix element of ^ J d 3 rh x j in a transversely po- 
larized state, [l3, 19]. This Fourier transform involves 
J\{Qb) and therefore the FRA corresponds to that of 
Eq. (Unj) and Eq. fl20|), with A = 1. Using this expan- 
sion, instead of simply taking the derivative of p 2l 
allows an expansion in basis functions that explicitly 



This paper is concerned with obtaining a general 
method to determine information about densities in 
the transverse plane. The use of Bessel series expan- 
sion, augmented by the finite radius approximation 
FRA of Eq. ©, Eq. ©, Eq. ^ and Eq. ^ al- 
lows us to determine the effects of experimental un- 
certainties and also allows us to estimate the effects 
of the incompleteness error caused by a lack of mea- 
surements at large values of Q 2 . The method can 
be applied to the extraction of any spatial quantity. 
One example, related to orbital angular momentum, 
is shown in Eq. and Eq. (|2H|) . 

The method is applied here to analyze electromag- 
netic form factors. We can see from Fig.[7]and Fig. [8] 
that the errors associated with the transverse charge 
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density and the two-dimensional Fourier transform 
of F2 are very small. The anomalous magnetiza- 
tion density pm, Fig. [SJ is also reasonably well de- 
termined, but future measurements extending our 
knowledge of F2 to higher values of Q 2 would reduce 
the existing uncertainties. 
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